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Verfahren zur rechnergestCitzten Ermittlung einer treibstoffoptimalen 

Ansteuerung von DUsen 

Die vorliegende Erfindung betrifft ein Verfahren zur rechnergestutzten Ermittlung 
einer treibstoffoptimalen Ansteuerung von Dusen. Ein solches Verfahren kann 
insbesondere zur Ansteuerung von Dusen eines Raumfahrzeuges wie belspielsweise 
eines Satelliten, einer Raumsonde, einer Raumstation oder Shnlichem Anwendung 
finden. 

Aus US 6,347,262 ist ein Verfahren zur treibstoffoptimalen Regelung von Lage und 
Nutation eines rotierenden Raumfahrzeuges bekannt. Hierbei wird Qber eine Bestim- 
mung von Drehgeschwindigkeiten und Drehwinkein ein Fehlersignals ermittelt und 
entsprechend dem Ergebnis fur dieses Fehlersignal ein Drehmoment auf das 
Raumfahrzeug ausgeubt. 

EP 0 977 687 B1 beschreibt unterschiedliche Moglichkeiten zur treibstoffoptimalen, 
rechnergestutzten Ansteuerung von Dusen eines Raumfahrzeuges. Es werden dabei 
uberwiegend Verfahren beschrieben, die einen Simplexalgorithmus beinhalten, 
wobei andererseits bereits dort ausgefuhrt wird, dass ein solcher Simplexalgorithmus 
mit einem hohen Aufwand an Rechenzeit verbunden ist. Das im Rahmen der 
' dortigen Erfindung beschriebene Verfahren benutzt ein Verfahren, das an den 
Simplexalgorithmus angelehnt ist und von der Bildung eines Simplextableaus 
ausgeht. Letztlich wird dort ein dualer Simplexalgorithmus zur Bildung eines 
treibstoffoptimalen Ansteuervektors verwendet. Als Alternativmoglichkeit zur Ermitt- 
lung einer treibstoffoptimalen Ansteuerung von DQsen wird in EP 0 977 687 B1 
lediglich ein ..Table Look-up"-Verfahren beschrieben. bei dem treibstoffoptimale 
DQsenanordnungen berechnet und in einer Tabelle abgelegt werden und aus diesen 
vorgespeicherten Ergebnissen fur den jeweils aktuellen Ansteuerfall ein aktuelles 
Ansteuerergebnis durch Kombination von gespeicherten Ergebnissen gebildet wird. 
Nachteilig ist jedoch, dass mit diesem Verfahren im Allgemeinen nicht eine 
tatsachlich treibstoffoptimale Losung gefunden wird. 
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Aufgabe der vorliegenden Erfindung ist es daher, ein Verfahren zur rechnergestutz- 
ten Ermittlung einer treibstoffoptimalen Ansteuerung von Dusen bereitzustellen. das 
mit einem moglichst geringen Rechenaufwand auskommt und trotzdem sicher zu 
einer treibstoffoptimalen Losung fuhrt. 

Diese Aufgabe wird gelost durch die Merkmale des Anspruchs 1. Weiterhin umfasst 
die Erfindung ein Computer-Progrannm nach Anspruch 5 und ein Computer- 
Programm-Produkt nach Anspruch 6. 

Die vorliegende Erfindung umfasst ein Verfahren zur rechnergestutzten Ermittlung 
einer treibstoffoptimalen Ansteuerung von DQsen gemSli einer Ansteuervorschrift 
b = Ax, wobei 

b einen gewunschten m-dimensionalen Krafte-ZDrehmomentvektor, 
A eine mxn-dimensionale Dusenmatrix und 

X den gesuchten n-dimensionalen Dusenansteuerungsvektor darstellt, mit 
Hilfe dessen ein entsprechender Schub durch die zugehorigen, angesteuerten Dusen 
bewirkt wird. Der Dusenansteuervektor soli dabei dem Minimierungskriterium 

genugen. Dieses Minimierungskriterium ist Ausdruck daftir, dass der Treibstoffver- 
braucli minimal sein soil, was durch eine Minimierung der Dusenansteuerung und 
damit durch eine Minimierung des durch die jeweiligen DQsen insgesamt erzeugten 
Schubes erzieit wird. 

GemdH der Erfindung sind nun Im Rahmen dieses Verfahrens folgende Verfahrens- 
schrltte vorgesehen: 

- rechnergestOtzt erfolgt eine definierte matrizentechnische Transformation von 
Ausgangsbedingungen fOr den Massenfluss der Dusen und des Minimierungs- 
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kriteriums, wobei die Ausgangsbedingungen und das Minimierungskriterium jeweils 
einer analogen matrizentechnischen Transformation unterzogen werden, 

- rechnergestutzt erfolgt eine datenverarbeitungstechnische Darstellung einer geo- 
metrischen Beschreibung der matrizentechnisch transformierten Ausgangsbedingun- 
gen, 

- durch ein rechnergestutztes, geometrisches Suchverfahren im Vektorraum erfolgt 
eine rechnergestutzte Ermittlung von Begrenzungs-Punktnnengen der geometrischen 
Beschreibung der Ausgangsbedingungen, und 

- das matrizenteclinisch transformierte Minimierungskriterium wird auf die Punkte der 
Begrenzungs-Punktmengen angewendet. 

Es kann also ein Simplexverfahren oder ein vergleichbares iteratives Verfahren im 
Rahman der Erfindung vermieden werden. Vielmehr wird statt eines iterativen 
Losungsweges ein rechnergestutzt realisierter geometrischer Losungsweg gewahlt. 
Durch die definierte matrizentechnische Transformation ist eine geometrische 
Beschreibung des Problems moglich, die ein geometrisches Auffinden einer treib- 
stoffoptimalen Losung eriaubt. Ein solches Verfahren kann viel schneller erfolgen als 
eine Losung des Problems uber ubiiche Simplexverfahren. 

Insbesondere kann im Rahmen der vortiegenden Erfindung vorgesehen sein, dass 

- zur matrizentechnischen Transformation der Ausgangsbedingungen fur den Mas- 
senfluss der Dusen eine homogene Losung der Ansteuervorschrifl gemali 

= 4,rdefiniert wird, wobei 

Ao\ die n x (n-m) dimensionale Nullraummatrix von A und 

n einen (n-m) dimensionalen Vektor beliebiger reeller Zahlen darstellt. 
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- im Rahmen der Anwendung der matrizentechnischen Transformation des Minimie- 
rungskriteriums eine rechhergestQtzte Berechnung von Skalarprodukten einer 
Vektordarstellung von Punkten der Begrenzungs-Punktmenge und des Vektors 



j=l 7=1 >1 



JP 



p:=n-m 



erfolgt und 



- eine treibstoffoptimale Losung mit Hilfe desjenigen Vektors r berechnet wird, 
dessen Skalarprodukt mit dem Vektor Vd minimal ist. 



Dies ist ein Beispiel, wie durch die Anwendung einer matrizentechnischen Trans- 
formation mit Hilfe der Nullraummatrix Ao sowohl das Minimerungskriterium wie auch 
die Ausgangsbedingungen geeignet transformlert werden konnen. Das Minimie- 
rungskriterium kann dabei im wesentlichen auf die einfache Bildung von 
Skalarprodukten von Vektoren, also auf eine einfache geometrische Rechenvor- 
schrifl, reduziert werden. 

Eine bevorzugte Weiterbildung der vorliegenden Erfindung sieht vor, dass 

- die matrizentechnisch transformierten Ausgangsbedingungen fQr den Massenfluss 
der Dusen rechnergestutzt in zulassige mehrdimensionale Wertebereiche umgerech- 
net werden, 



- zur Ermittlung der Begrenzungs-Punktmengen eine Bildung von mindestens einer 
mehrdimensionalen Schnittmenge der einzelnen zulassigen mehrdimensionalen 
Wertebereiche erfolgt und 

- die Begrenzungs-Punktmengen als diejenigen Punktmengen ermittelt werden, die 
die mindestens eine Schnittmenge begrenzen, 

Auf diese Weise kann die Anzahl der im Rahmen des erfindungsgemalien Verfah- 
rens zu betrachtenden Punkte deutlich reduziert werden, da lediglich diejenigen 
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Punkte betrachtet werden, die die mindestens eine Sclinittmenge begrenzen, Diese 
bedeutet eine weitere Verkurzung der fur das Verfahren benotigte Rechenzeit und 
damit einen weiteren Vorteil gegenuber den bisher bekannten Verfahren. 

Das vorgenannte Verfahren kann vorteilhaft dadurch noch weiter vereinfacht und 
damit die erforderiiche Rechenzeit weiter reduziert werden, dass 

- eine wiederholte Projektion der zuiassigen nnehrdimensionalen Wertebereiche der 
Dimension p auf eine Dimension p-1 erfolgt, bis eine Projektion der zulSssigen 
Wertebereiche auf Begrenzungsintervalle einer Dimension p=1 erreicht ist und 

- anschlieRend ein rechnergestutztes Suchverfahren eine rechnergestutzte Ermitt- 
lung von Begrenzungs-Punktmengen als Schnittmenge von Begrenzungsintervallen 
durchfuhrt. 

Durch die hier vorgesehene wiederholte Projektion von einer Dimension p auf eine 
Dimension p-1 gelangt man also schrittweise zu immer niedrigeren Dimensionen (p- 
1, p-2, p-3 usw.), in denen die fraglichen Wertebereiche dargestellt werden, und 
damit zu rechentechnisch einfacher handhabbaren Problemen, als sie es in der 
Ausgangsdimension p waren. Grundsatzlich konnte bei Erreichen einer geeigneten, 
ausreichend reduzierten Dimension p>Pr>1 bereits die Ermittluhg der Begrenzungs- 
Punktmengen erfolgen. Bevorzugt wird jedoch die wiederholte Projektion so lange 
durchgefOhrt, bis eine Beschreibung der Dimension p=1 erreicht ist. Hier kann das 
Problem auf einfachste Weise und mit dem geringsten Aufwand an Rechenzeit 
gelost werden. 

Ein weiterer Gegenstand der vorliegenden Erfindung ist ein Computer-Programm zur 
rechnergestutzten Ermittlung einer treibstoffoptimalen Ansteuerung von Dusen 
gemali einer Ansteuervorschrifl b = Ax, wobei 

b einen gewunschten m-dimensionalen Krafte-ZDrehmomentvektor, 
A eine mxn-dimensionale DUsenmatrix und 

X den gesuchten n-dimensionalen Dusenansteuerungsvektor darstellt, und 
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der Dusenansteuervektor dem Minimierungskriterium 
geniigen soil. 

GemaB der Erfindung ist vorgesehen, dass das Computerprogramm folgendes 
beinhaltet: 

- eine erste Programmroutine zur rechnergestutzten Durchfuhaing einer definierten 
matrizentechnischne Transformation von Ausgangsbedingungen fur den Massen- 
fluss der Dusen und des MinimierungskriteriumSi 

- eine zweite Programmroutine zur rechnergestutzten Durchfuhrung einer datenver- 
arbeitungstechnischen Darstellung einer geometrischen Beschreibung der matrizen- 
technisch transformierten Ausgangsbedingungen. 

- eine dritte Programmroutine zur rechnergestutzten Durchfuhrung eines geometri- 
schen Suchverfahrens im Vektorraum zur rechnergestutzten Ermittlung von Begren- 
zungs-Punktmengen der geometrischen Beschreibung der Ausgangsbedingungen, 

- eine vierte Programmroutine zur rechnergestutzten Anwendunig des matrizen- 
technisch transformierten l\4inimierungskriteriums auf die Punkte der Begrenzungs- 
Punktmengen, 

Ein solches Computer-Programm ist zur Durchfuhrung des vorgenannten erfindungs- 
gemafien Verfahrens geeignet. Es konnen auch weitere Programmroutinen im 
Rahmen dieses Computer-Programms vorgesehen werden, die zur Durchfuhrung 
einer oder mehrerer der vorgenannten Weiterbildungen des erfindungsgemaBen 
Verfahrens geeignet sind. 

Ein weiterer Gegenstand der vorliegenden Erfindung ist ein Computer-Programm- 
Produkt, beinhaltend einen maschinenlesbaren Programmtrager, auf dem ein vorste- 
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hend beschriebenes Computer-Programm in Form von elektronisch auslesbaren 
Steuersignalen gespeichert ist. Die Steuerslgnale kdnnen in jeder geeigneten Form 
gespeichert sein, die elektronische Auslesung kann dann entsprecliend durch 
elektrische, magnetische, elektromagnetische, elektrooptische oder sonstige elektro- 
nische Verfahren erfolgen. Beispiele fur solclie Programmtrager sind Magnetbander, 
Disketten, Festplatten, CD-ROM oder Halbleiterbausteine. 

Ein spezielles Ausfuhmngsbeispiel der vorliegenden Erfindung wird nachfolgend 
aniiand der Figuren 1 bis besclirieben. 

Es zeigen 

Fig. 1 : Zulassiger Bereicli fQr die Vektoren r in einem Raum der Dimension p=2, 
wobei eine Bedingungsungleicliung rii^ir-Z^n.) > 0 erfullt ist 

Fig. 2: Zulassiger Bereicli fur r im zweidimensionalen Raum unter Beachtung der 
Ausgangsbedingungen fQr mehrere DGsen 

Fig. 3: Zulassiger Bereich fur r im eindimensionalen Raum unter Beacfitung der 
Ausgangsbedingungen fur mehrere DQsen 

Fig. 4: Schnittmenge zweier (p-1 )-dimensionaler Ebenen 

Fig. 5: Lineare Koordinatentransformation von Einheitsvektoren 

Die Erfindung betrifft das folgende Problem: Es sollen durch Dusen eines 
Raumfahrzeuges gewunschte Krafte und/oder Drehmomente auf das Raumfahrzeug 
ausgeubt werden. Wie miissen die Dusen angesteuert werden, damit die 
gewunschten Krafte und/oder Drehmomente mit einer minimalen Menge an 
Treibstoff fur die Dusen erzielt werden konnen? Zusatzlich muss der Schub der 
DQsen in einem eriaubten Wertebereich liegen, also zwischen einem minimal 
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moglichen und einem maximal mOglichen Wert. Entsprechend mussen daher auch 
die Ansteuerwerte sich in einem eriaubten Wertebereich bewegen. 

Dieses Problem ist ein sogenanntes lineares Optimierungsproblem. Bislang wurden 
zur Losung solcher Probleme meist Simplex-Algorithmen verwendet, wie eingangs 
beschrieben, die aber die ebenfalls eingangs gescinilderten Nachteile aufweisen. Die 
vorliegende Erfindung stelit eine verbesserte Moglichkelt zur L6sung des linearen 
Optimierungsproblems dar. 

Das lineare Optimierungsproblem fur eine treibstoffoptimale Dusenansteuerung 
beruht auf der Gleichung: 

b=^Ax (2.1) 

wobei 

b einen gewunschten m-dimensionalen Krafte-ZDrehmomentvektor, 
A eine mxn-dimensionale Dusenmatrix und 

X den gesuchten n-dimensionalen Dusenansteuerungsvektor darstellt. 

Es soil fur den gesuchten, unbekannten Vektor x ein Minimierungskriterium erfQIlt 
sein, durch das eine treibstoffoptimale Ansteuerung garantiert wird, nSmlich 

(2.2) 

Weiterhin mussen sich die Werte fiir den Ve^or x in einem zulassigen Wertebereich 
bewegen analog dem zulassigen Wertebereich fur den Schub der Dusen: 

0<Xi <lfori = l,...7i . (2,3) 

Ohne BeschrSnkung der Allgemeinheit des Problems dargestellt in (2.1) kann 
angenommen werden, dass die Anzahl der Dusen n grdlier ist als die Anzahl m der 
gewunschten Krafte und/oder Drehmomente 

n>m (2.4) 
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und dass die Dusenmatrix >4 vollen Rang hat, 

rank(A)=m 



(2.5) 



Ein VerstoR gegen (2.5) wOrde beispielsweise bedeuten, wenn alle Dusen exakt in 
die gleiche Richtung weisen wurden. Dies wird jedoch stets durch eine sinnvolie 
5 Anordnung der Dusen an dem Satelliten vermieden. 

Es kann nun das Minimierungskriterium umgeschrieben warden: 
.Die allgemeine Losung der Gleichung (2.1) kann wie folgt geschrleben werden 

^g=Xpa-^XHo (2.6) 



Eine spezielle Losung Xpa kann folgendermaBen aus (2.1) generiert werden: 
15 b = Ax=^ AA'iAA^y'Ax^ AA'iAA'y'b = AiA'iAA'y'b) = A^x^) (2.7) 

Daraus erhalt man die folgende spezielle Losung: 



( 20 Die allgemeine homogene Losung des Problems bilden alle Vektoren Xho fur die gilt: 
^ ^*o=0, (2.9) 

Diese homogene Losung kann folgendermaBen geschrleben werden: 




wobel 



spezielle Losung von (2.1) 
Xho allgemeine Losung der homogenen Gleichung von (2.1 ). 



(2.8) 



(2.10) 



25 



wobel 



Ao: 



n X (n-m) Nullraummatrix von A 

(n-m) dimensionaler Vektor beliebiger reeller Zahlen. 



n 
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Das Minimierungskriterium aus (2,2) kann nun mit Hilfe der Gleichungen 
(2.6,2.8,2.10) geschrieben werden als 

•/W = l;(^pa)y+v/r (2.11) 

5 mit 

(2.12) p:=n-m. 

(2.13) 



:= 



7=1 7=1 y=l 



Gleichung (2.11) kann nun folgendemiafien interpretiert werden: Fur die eriaubten 
10 Werte von rsoll derjeinige Vektor gefunden werden, dessen Skalarprodukt mit dem 
Vektor minimal wird. Dies folgt aus der Tatsache, dass die spezielle Losung wie in 
Gleichung (2.8) definiert einer weiteren Minimierung nicht mehr zuganglich ist. Die 
vorliegende Erfindung setzt diese Suche nach dem opimalen Vektor r 
rechnergestutzt mit Hilfe einer entsprechenden Programmroutine eines Gomputer- 
15 Programms um. 



c > 



Geometrische Beschreibung der treibstoffoptimalen Losung: 

Die zulassigen Werte bzw. der zulassige Bereich fur r wird aus den 
Ausgangsbedingungen fur den Massenfluss der Dusen nach (2.3) durch eine 
'^20 matrizentechtnlsche Transformation gewonnen. Man erhSIt dies formal durch 
Einsetzen von (2.6,2.8,2.10) in 2.3) 

-x^<^r<l-x^, (2.14) 

wobei 

25 1: n-dimensionaler Vektor, bei dem alle Elemente gleich 1 

sind 

Es kann nun eine geometrische Beschreibung der so transformierten Ausgangsbe- 
dingungen erfolgen, die im.Rahmen der vorliegenden Erfindung datenverarbeitungs- 
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technisch realisiert wird. Es kann dabei die Gleichung (2.14) mit den folgenden 
Definitionen in skalarer Fomn geschrieben werden. 



4.= 



*01 



*02 



a, 



On 



Damit erhalt man: 



fori=1,..n 
fori=1,..n. 



(2.15a-c) 



(2.16) 
(2.17) 



DIvidiert man nun (2.16, 2.17) durch \ao,\ - aufgrund der Bedlngung (2.5) ist dies 
immer moglicli - und mit den folgenden Definitionen 



10 



15 



n.. := 



I ^0/ I 



fori=1,...n 



On 



n. := 5£_ for i=n+1,...2n 

I «0/ I 

A, := fori=1,...n 

I «0/ I 



A, :=• 



ioo/ I 



fori=n+1,...2n 



(2.18) 
(2.19) 
(2.20) 
(2.21) 



konnen die Gleichungen (2.16, 2.17) wie folgt geschrieben werden: 

/i,''(r->l../i/)^0 fori=1....2n. (2.22) 



Damit kann nun fOr einen bestimmten Wert von / die Gleicliung (2.22) als eine ein- 
20 Oder melirdimensionale Ebene interpretlert und damit entsprechend datenverarbel- 
tungsteciiniscii dargestellt werden (sielie Fig. 1), wobei diese Ebene um den Vektor 
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Xn, bezQglich des Ursprungs verschoben ist und senkrecht zum Nomialenvektor n, 
orientiert ist. Im Fall der Fig. 1 ist diese ein- oder mehrdimenslonale Ebene aus 
Grunden der Einfachheit als Gerade dargestellt, also als eine Ebene der Dimension 
1 . Diese eindimensionale Ebene begrenzt den Bereich der zulSssigen Werte bzw. 
der zulassigen Vektoren fOr r. Wenn das Vorzelten von negativ ist, ist auch der 
Nullvektor in der eriaubten Region fQr r enthalten, ist jedoch das Vorzeichen von Xj 
positiv, dann ist der Nullvektor aus der eriaubten Region ausgeschlossen. Der 
' Normalenvektor n, zeigt in Richtung der eriaubten Region. Die geometrische 
Beschreibung der Gleichung (2.22) ist fQr p=2 in Fig. 1 fQr Xi>0 dargestellt. Es wird 
I darauf hingewiesen, dass zu jeder Bedingung / in (2.22) eine Gegenbedingung i+n 
mit n;+n=-/7/ existiert (entsprechend den Ausgangsbedlngungen fQr den minmalen und 
maximalen Massenfluss einer individuellen Duse), woraus n Paare von Bedingungen 
entstehen die zulassige „Gebietsstreifen" definieren, wie in Fig. 2 dargestellt. 

Die gemeinsame Schnittmenge aller zulassigen Bereiche bestimmt denjenigen Be- 
reich im p-dimensionalen Raum, in dem r alle 2n transformierten Ausgangsbedln- 
gungen bzw. alle Bedingungen aus (2.22) erfullt. In Fig. 2 ist dies fur n=3 DQsen, also 
fur 6 Bedingungen, dargestellt. Wenn die Schnittmenge aus alien zulassigen 
Bereichen gebildet wird, ergibt sich ein konvexer Bereich als Schnittmenge, d.h. dass 
fur jeweils zwel Punkte, die innerhalb der Schnittmenge liegen, auch alle Punkte auf 
einer geraden Linie zwischen diesen beiden Punkten innerhalb der Schnittmenge 
liegen. 

Aber nicht nur der Bereich der zulassigen Werte fQr r kann geometrisch beschrieben 
werden, sondem auch der optimale Wert fur r, ropt, welcher J in Gleichung (2.11) 
minimiert: J wird, wie gesagt, minimal fur denjenigen Vektor r der das kleinste 
Skalarprodukt mit oder die kleinste Projektion auf den Vektor Vd aufwelst. In Fig. 2 ist 
dies fur den Vektor r gegeben, der zu dem mit einem Kreis markierten Eckpunkt, 
welcher auf der Geraden g, liegt, zu der der Nomialenvektor ni gehdrt. Der optimale 
Punkt ergibt sich also als einer der Eckpunkte der Begrenzungslinien der 
Schnittmenge. Die Begrenzungslinien bilden also Begrenzungs-Punktmengen. 
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Die geometrische Beschreibung dieses zweidimensionalen Beispiels kann in 
analoger Weise unmittelbar auf hohere Dimensionen p angewendet werden, 

Wenn das Problem ein dreidimensionales Problem ist, dann konnen die n transfor- 
5 mierten Ausgangsbedinungen als n Gebietsbereiche zwischen n Ebenenpaare 
beschrieben werden, wobei jedes Ebenenpaar aus zwei zueinander parallelen 
Ebenen besteht und der dazwischenliegende Gebietsbereich den Bereich fQr r 
darsteilt, der die zugehorigen transformierten Ausgangsbedingungen erfullt. Der 
/-insgesamt zulassige Bereich fur r, der alie transformierten Ausgangsbedingungen 
C erfullt, ergibt sich als Schnittmenge aller n Gebietsbereiche und entspricht damit 

f ^ ' W einem dreidimensionalen Polygon. Der optimale Punkt ergibt sich auch hier als einer 
der Eckpunkte der Begrenzung des Polygons, welches nun durch Begrenzungs- 
flachen als Begrenzungs-Punktmengen begrenzt wird. Der optimale Punkt ist 
wiederum derjenige Eckpunkt, dessen Vektor r mit dem Vektor Vd das kleinste 
15 Skalarprodukt aufweist. 

Wenn es sich urn ein eindimensionales Problem handelt, dann konnen die Bereiche, 
die durch die n transformierten Ausgangsbedingungen begrenzt werden, als 
eindimensionale Intervalle beschrieben werden. Die n transformierten Ausgangsbe- 
20 dingungen werden dann durch Begrenzungspunkte beschrieben. Dies ist beispielhaft 
( i^^J^ '^"9- 3 dargestellt Die in Fig. 3 dargestellten Klammern gleicher GroSe bilden 



Wrfn'eweils ein Paar von transformierten Ausgangsbedingungen. Die unteren Grenzen 
der Ausgangsbedingungen sind mit ^3. bezeichnet, die oberen Grenzen der 
Ausgangsbedingungen mit U,... (Gleichungen 2.20, 2.21). Im Beispiel der Fig. 3 
25 existiert ein Schnittmengen-lntervall als zulassiger Bereich fur r, das gleichzeitig alle 
transformierten Ausgangsbedingungen erfQIIt. Der optimale Wert fur r ergibt sich 
wiederum aus einem Punkt aus der Begrenzungs-Puntkmenge, hier als der rechte 
Begrenzungspunkt des Schnittmengen-lntervalls. Der zu diesem Punkt gehorende 
Vektor r hat wiederum das kleinste Skalarprodukt mit dem ebenfalls in Fig. 3 
30 dargestellten, negativen Vektor Vd. 
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Die gleichen Interpretationen und Beschreibungen sind auf Dimensionen p>3 
anwendbar. 

Rechnergestutzte Ermittlung einer treibstoffoptimalen Losung: 

Es soli nun ein spezielles Beispiel zur Ermittlung einer treibstoffoptimalen Losung 
dargestellt werden. Es wird dabei ein Verfahren wie im Rahmen der Erfindung 
beschrieben durchgefuhrt, d. h. es erfolgt insbesondere rechnergestOtzt eine 
Ermittlung von (p-l)-djmenslonalen Begrenzungs-Punktmengen als geometrische 
Beschreibung der transformierten Ausgangsbedingungen in Form von Ebenen. 
Dieses spezielle Beispiel beschreibt nun eine bevorzugte Weiterbildung der 
Erfindung und beinhaltet Insbesondere die folgenden Verfahrensscfiritte: 

1. RechnergestQtzte Bestimmung einer ersten der ermittelten Ebenen als 
Begrenzungs-Punktmenge, die wahrschelnlich die treibstoffoptimale Losung 
beinhaltet. Weltere Detail zum Auffinden dieser Ebene siehe weiter unten. 

2. Rechnergestutzte Anwendung einer linearen Koordinatentransformation auf 
alle Vektoren r so dass der Normalenvektor der ersten Ebene anschlieliend 
parallel zu dem p-ten Einheltsvektor derjenlgen Einheitsvektoren ist, die den 
Vektorraum der Vektoren r aufspannen, und so dass der Normalenvektor in 
die gleiche Richtung weist wie dieser Einheltsvektor. 

3. RechnergestQtzte Berechnung von Schnittmengen in Form von Schnlttebenen 
der ubrigen ennlttelten Ebenen mit der ersten Ebene. Die p-te Komponente 
jeder Schnlttebene hat dabei denselben konstanten Wert wie die p-te 
Komponente der ersten Ebene nach der linearen Koordinatentransformation. 
Damit kann die eindeutige geometrische Beschreibung der Schnlttebenen 
vereinfacht werden auf eine Dimension Pnew= p-l- Damit sind zur 
Beschreibung eines Vectors Vd nur noch die ersten p-l Komponenten 
erforderlich. 
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4, Rechnergestutzte PrQfung: ist pnew verschieden von 1? Dann werden die 
Schritte 1-3 emeut durchgefQhrt. 

RechnergestQtzte Priifung: ist Pnew=1 dann erfolgt eine rechnergestutzte 
Prufung, ob ein Intervail als Schnittmenge vorliegt, das einen zulassigen 
Bereich fur rdarstellt. Wenn kein solches Intervail vorliegt, wird das Verfahren 
beendet, gemafi Schritt 1 eine zweite Ebene ausgewahit und das Verfahren 
emeut durchgefQhrt. Liegt ein solches Intervail vor, so erfolgt eine 
rechnergestutze Multiplikation der Begrenzungspunkte des Intervalls mit dem 
mittlerweile zu einem Vektor der Dimension 1, also zu einem Skalar 
reduzierten Vektor Vd. Es wird der kleinere Wert der Ergebnisse der 
Multiplikation abgespeichert, und falls zuvor bereits das Verfahren fur eine 
Oder mehrere andere Ebenen durchgefQhrt wurde, mit zuvor gespeicherten 
Werten verglichen. Wenn der zuletzt gespeicherte Wert kleiner Ist als ein 
zuvor gespelcherter Wert, dann wird das Verfahren erneut gemSR Schritt 1 mit 
einer neuen Ebene begonnen und entsprechend durchgefQhrt. Wenn der 
zuletzt gespeicherte Wert nicht kleiner ist, dann reprasentiert der zuvor 
gespeicherte Wert die optimale Losung. 

Es sollen nun die einzelnen rechnergestQtzten Schritte detaillierter betrachtet 
werden. 

Bestimmung der ersten der ermittelten Ebenen nach Schritt 1 : 

Grunds3tzlich kann die Ebene, die die treibstoffoptimale Losung beinhaltet, nicht a 
priori eindeutig ermittelt werden. Am Beispiel des Falles p=2 kann jedoch verdeutlicht 
werden, dass bestimmte grundsatzliche Aussagen uber die gesuchte (p-1)- 
dimensionale Ebene bzw. Gerade (im Fall p=2 wird besitzt also die gesuchte Ebene 
die Dimension 1) mdglich sind. Aus geometrischen Oberlegungen wird klar, dass die 
treibstoffoptimale Losung immer auf einem der Geradenabschnitte liegt, welche die 
Schnittmenge aller zulassigen Bereiche begrenzen, und zwar auf derjenigen 
Geraden, deren Normalenvektor (der in Richtung des zulassigen Bereiches weist) 
das grolite Skalarprodukt mit dem Vektor Vd besitzt und dessen Richtung damit am 
besten der Richtung des Vektor Vd entspricht. In Fig. 2 entspricht der Normalenvektor 
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ni am besten dem Vektor va so dass der Vektor Vd mit diesem Normalenvektor das 
grolite Skalarprodukt aufweist. Dieses Skalarprodukt kann rechnergestutzt berechnet 
werden und der entsprechende Normalenvektor und damit die zugeliorige (p-1)- 
dimensionale Ebene rechnergestutzt bestimmt werden. Bei diesem Verfahren wird in 
Bezug auf Fig. 2 also davon ausgegangen, dass die treibstoffoptimale Losung auf 
dem Geradenabsclinitt liegt, der die Sclinittmenge begrenzt und zu des Gerade 
durch den Normalenvektor ni charaktierisiert wird. Analog wird in alien Anwendungs- 
fallen dieses Verfahrens vorgegangen. 

Es werden also insbesondere die folgenden Verfahrensschrltte durchgefuhrt: 

• Bestimmung desjenlgen Normalenvektors, dessen Skalarprodukt mit dem 
Vektor Vd maximal ist. 

Ermittlung, ob die zugehorige Gerade (oder allgemein: (p-l)-dimensionale 
Ebene) eine Begrenzungs-Punktmenge der Schnlttmenge darstellt (s. unten). 
Falls nein, so wird diese Gerade verworfen und diejenige Gerade herange- 
zogen, deren Normalenvektor das zweitgrofite Skalarprodukt mit dem Vektor 
Vd aufweist. Dann wird mit dieser Gerade erneut gepruft, ob sie eine 
Begrenzungs-Punktmenge darstellt. 

Koordinatentransformation der ersten Ebene: 

Urn die rechnergestutzte Ermittlung der Schnittebenen zu vereinfachen, welche sich 
als Schnittmenge der einzelnen (p-1 )-dimensionalen Ebenen ergibt, wird ein neues 
Koordinatensystem eingefuhrt. Dies 1st in Fig. 4 dargestellt. Das Koordinatensystem 
wird durch eine lineare Koordinatentransformation so gedreht, dass ein Einheits- 
vektor (es wird hier derjenige Einheitsvektor mit dem hochsten Index, also mit Index 
p gewahit) mit dem Normalenvektor derjenigen Ebene zusammenfailt, die im 
vorstehend beschriebenen Schritt als erste Ebene ausgewahit wurde. Als Vorteil 
ergibt sich durch diese Koordinatentransformation, dass die p-te Koordinate aller 
Schnittebenen nun den gleichen konstanten Wert aufweisen. In Fig. 4 sind zwei 
Ebenen pi und p2 dargestellt, die im Fall der Fig. 4 eine Schnittebene gi2 der 
Dimension 1, also eine Gerade als Schnittmenge aufweisen. Es sei nun pi die nach 
dem vorstehend beschriebenen Verfahren ausgewahlte erste Ebene. Dann wird 
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durch eine Koordinatentransformation der Einheitsvektor 63 urn 180° urn den Vector 
©3 + rti gedreht, urn einen transformierten Einheitsvektor 63 zu erhaiten, der nun mit 
dem Vektor ni zusammenfallt (siefie Fig. 5). Die dritte Koordlnate von gr2 wird dann 
eine Konstante. Diese Koordinatentransformation kann ausgedriickt werden als: 
C3' = n, = [idd"^ - E^, (2.23) 

mit 



(2.24) 



Die Transformationsgleichung (2.23) kann von drei Dimensionen auf ein p- 
dimensionales Problem verallgemeinert werden mit: 

e^' =n,=[2dd^ -E^^ (2.25) 

mit 

e„ +n, 



d:= 



wobei / der Index der gewahlten Ebene ist. 



(2.26) 



Bestimmung der Schnittebenen / Schnittmengen: 

Es konnen Gleiciiungen zur Beschreibung der Schnittmengen in Form von Ebenen- 
gleichungen ermittelt werden. Die Ebenengleichung fQr jede Ebene / im p-dimensio- 
nalen Raum kann geschrieben werden als: 



'p-i 



-A; =0 



(2.27) 



Die Schnittmenge dieser Ebene mit einer gewShlten Ebene ; wird beschrieben durch 
die Gleichung: 

(2.28) 
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Damit erhalt man als Beschreibung der Schnittmenge in Fomn einer Schnittebene: 



+ 71,^ Ay - A, =0. 



(2.29) 



Gleichung (2.29) kann in Form der Gleichung (2.27) umgesclirieben werden zu 



(2.30) 



mit 



«, := 



n 



n 



n 



ip-i 



r := 



(2.31) 



(2.32) 



(2.33) 



Rechnergestutztes Suchverfaliren: 

Das rechnergestutzte Suchverfahren kann nun grundsdtzlich folgende Schritte 
beinhalten: 



1. Rechnergestutzte Ermittlung der erfolgversprechendsten ersten Ebene als 
Startebene, welcher die treibstofFoptimale Losung am wahrscheinlichsten 
beinhaltet und Speicherung des Vektors, der vom Koordinatenursprung zu 
dieser Ebene weist. 

2. Rechnergestutzte Durchfuhrung einer Koordinatentransfonmation wie oben 
beschrieben und rechnergestutzte Ermittlung von Schnittebenen, wie 
ebenfalls oben beschrieben. 

3. Mit den ermittelten Schnittebenen werden dann die vorstehend beschriebenen 
Verfahrensschritte durchgefUhrt, bis die Dimension auf p=1 reduziert ist. 
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4. Es wird dann ermittelt ob ein zulassiges Intervall vorliegt (siehe oben). 1st dies 
gegeben, wird entsprechend der obigen Beschreibung ein Punkt als 
Beschreibung einer treibstoffoptimalen Losung ermitttelt und durch sine 
entsprechende RQclctransformation der Koordinaten die treibstoffoptimale 
Losung rechnergestutzt berechnet. 

5. Wiederholung der Scliritte 1 bis 4, bis keine weitere Verbesserung mehr 
erzielt wird. Die dann vorliegende LOsung stellt die globaie treibstoffoptimale 
Losung dar. 

Fig. 6 zeigt inn Ralimen eines Berechnungsbeispiels einen Vergleicli zwischen der 
Rechenzeit eines ubiichen Verfahrens M1 nach einem Sinnplex-Algorithmus mit dem 
erfindungsgennafien Verfahren M2, Jewells fiir Anordnungen mit 4, 5, 6, 7 oder 8 
Dusen. Wie bereits aus Fig. 6 a) deutiicfi wird, liegt in diesem Beispiel die Rechenzeit 
fiir das erfindungsgemaBe Verfahren M2 fur die Falle mit 4, 5 und 6 Dusen deutlich 
unter der Rechenzeit des Simplex-Verfahrens M1. Die Fig. 6 b) veranschaulicht das 
Verhaltnis Rechenzeit M1 / Rechenzeit M2. Fur 4 Dusen ist die Rechenzeit von M2 
um einen Faktor 100, fur 5 Dusen um einen Faktor 16 und fQr 6 Dusen urn einen 
Faktor 3 langsamer als das erfindungsgemSBe Verfahren. Erst fur hohere 
Dusenzahlen wird in diesem Beispiel das erfindungsgemSlie Verfahren M2 
langsamer als das Simplex-Verfahren. 

Soli die optimale LOsung im eindimensionalen Raum ermittelt werden (d.h. ist die 
Differenz zwischen der Zahl der Dusen und der Zahl der Bedingungen gleich 1) so 
kann die oben ausgefuhrte Darstellung des Verfahrens in eine kompaktere Form 
gebracht werden. Dann ergibt sich Gleichung (2.14) als n Gleichungen fur die 
Obergrenze und Untergrenze, also den zulassigen Bereich der transformierten 
Ausgangsbedingungen 

-x^<aor<l-x^ (2.23) 

mit dem n-dimensionalen Nullraumvektor ao und dem unbekannten Skalar r. 

Wenn eine Kompontente ao,- gleich Null ist, dann muss fur die spezielle Ldsung 

0<=Xpa/<=1 gelten, ansonsten kann die o.g. Bedingung nicht erfQIlt werden. Dieser 
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Fall ist jedoch hier ausgeschlossen, da dieser Fall in der Praxis eine schlecht 
gewahlte Dusenanordnung bedeuten wQrde. 

Die Gleichungen (2.23) werden durch die Komponenten ao/ (unter Beachtung des 
Vorzeichens von aoi) dividiert, und es ergibt sich 



mit 



''mm 


max 








r ' 

mm/ 


















r 

maxi 







(2.24) 



fora^i >0 
fora^, <0 

fora^t >0 
fora^. <0 



(2.25) 



(2.26) 



Die Bedingungen in (2.24) kOnnen nur dann gleichzeitig erfOllt werden, wenn 

max(r^ ) < min(r^ ) (2.27) 
Die optimaie Losung fur r , ropt, hangt von dem Vorzeichen des Skalars Vd in dem 
Minimierungskriterium (2.11) ab: Fur va grolier Null, dann ist das 
iVIinimierungskriterium so klein als moglich, wenn r am linken Rand des Intervalls 
(2.27) gewahit wird. 

r^, = max(r^ )fiirvj>0 (2 .28) 

und andererseits gilt fUr Vd kleiner Null, dass das Minimierungskriterium so klein als 
mOglich , wenn r am rechten Rand des Intervalls (2.27) gewahit wird, 

rop, = min(r^ )Jurvj<0 (2.29) 

Fur \/d=0 kann jeder beliebige Wert fur r aus dem Interval (2.27) gewahit werden, um 
das Minimierungskriterium zu erfullen. 

Die Werte fur ao und va hangen nur von der Position und der Schubrichtung der 
DQsen ab, aber nicht von den aktuellen Schubanforderungen. Daher muss die 
Fallunterscheidung gemaii (2.25, 2.26. 2.28, 2.29) nur ein einziges Mai vor der 



TU/bg; 13.3.2003 



•P610954/DE/1 ' /' 21 



rechnergestdtzten Durchfuhrung des Verfahrens zur Ermittlung der optimalen LOsung 
erfolgen. 
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Patentanspruche: 

1. Verfahren zur rechnergestutzten Ermittlung einer treibstoffoptimalen Ansteue- 
rung von Dusen gemafi einer Ansteuervorschrift Ax, wobei 
b einen gewunschten m-dimensionalen Krafte-ZDrehmomentvektor, 
A eine mxn-dimensionale Dusenmatrix und 

X den gesuchten n-dimensionalen Dusehansteuerungsvektor darstellt, und 
der Dusenansteuervektor dam Minimierungskriterium 

^min 

genugen soli, 

dadurch gekennzeichnet, dass 

- rechnergestutzt eine definierte matrizenteciinische Transformation von Aus- 
gangsbedingungen fur den Massenfluss der Dusen und des Minimierungs- 
kriteriums erfolgt, 

- rechnergestQtzt eine datenverarbeitungstechnische Darstellung einer geo- 
metrischen Beschreibung der matrizentechnisch transformierten Ausgangs- 
bedingungen erfolgt, 

- durch ein rechnergestutztes, geometrisches Suchverfahren inn Vektorraum 
eine rechnergestutzte Ermittlung von Begrenzungs-Punktmengen der 
geometrischen Beschreibung der Ausgangsbedingungen erfolgt, und 

- das matrizentechnisch transformierte Minimierungskriterium auf die Punkte 
der Begrenzungs-Punktmengen angewendet wird. 

2. Verfahren nach Anspruch 1 , 
dadurch gekennzeichnet, dass 

- zur matrizentechnischen Transformation der Ausgangsbedingungen fur den 
Massenfluss der Dusen eine homogene Losung der Ansteuervorschrift 
gema& = yi^r definiert wird, wobei 

Aq: die n x (n-m) dimensionale Nullraummatrix von A und 

n einen (n-m) dimensionalen Vektor beliebiger reeller Zahlen darstellt, 

- im Rahmen der Anwendung der matrizentechnischen Transformation des 
Minimierungskriteriums eine rechnergestutzte Berechnung von Skalarpro- 
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dukten einer Vektordarstellung von Punkten der Begrenzungs-Punktmenge 
und des Vektors 



:= 



p'j=n — m 



.7=1 i=i 

erfplgt und 

eine treibstoffoptimale Losung mit Hilfe desjenigen Vektors r berechnet wird, 
dessen Skalarprodukt mit dem Vektor v^i minimal ist. 



3. Verfahren nach einem der Anspruche 1 Oder 2, 
dadurcli gekennzeichnet, dass 

- die matrizenteclinisch transformierten Ausgangsbedingungen fur den 
Massenfluss der Dusen rechnergestUtzt in zulSssige mehrdimensionale 
Wertebereiche umgerechnet werden, 

- zur Ermittlung der Begrenzungs-Punktmengen eine Bildung von mindestens 
einer melirdimensionalen Schnittmenge der einzelnen zulassigen mehrdi- 
mensionalen Wertebereiclie erfolgt und 

- die Begrenzungs-Punktmengen als diejenigen Punktmengen ermittelt wer- 
den, die die mindestens eine Schnittmenge begrenzen, 

4. Verfahren nach Anspruch 3, 
dadurch gekennzeichnet, dass 

- eine wiederholte Projektion der zulassigen mehrdimensionalen Wertebe- 
reiche der Dimension p auf eine Dimension p-1 erfolgt, bis eine Projektion 
der zulassigen Wertebereiche auf Begrepzungsintervalle einer Dimension 
p=1 erreicht ist und 

- anschlieUend ein rechnergestutztes Suchverfahren eine rechnergestutzte 
Ermittlung von Begrenzungs-Punktmengen als Schnittmenge von Begren- 
zungsintervallen durchfuhrt. 



5. Computer-Programm zur rechnergesttitzten Ermittlung einer treibstoffoptima- 
len Ansteuerung von Dusen gemali einer Ansteuervorschrift b = Ax, wobei 
b einen gewunschten m-dimensionalen.Krafle-/prehmomentvektor, 
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A eine mxn-dimensionale Dusenmatrix und 

X den gesuchten n-dimensionalen DQsenansteuerungsvektor darstellt, und 
der Dusenansteuervektor denri Minimierungskriterium 

J := 2!C^i ->min 

5 genugen soil, 

dadurch gekennzelchnet, dass das Computerprogramm beinhaltet 

- eine erste Programmroutine zur rechnergestQtzten Durchfuhrung einer 
definierten matrlzentechnischne Transformation von Ausgangsbedingungen 
fur den Massenfluss der DQsen und des Minimierungskriteriums, 

il - eine zweite Programmroutine zur rechnergestQtzten Durchfuhrung einer 

datenverarbeitungstechnischen Darstellung einer geometrischen Beschrei- 
bung der matrizentechnisch transformierten Ausgangsbedingungen, 

- eine dritte Programmroutine zur rechnergestQtzten Durchfuhrung eines 
geometrischen Suchverfahrens im Vektorraum zur rechnergestQtzten 

15 Ermittlung von Begrenzungs-Punktmengen der geometrischen Beschreibung 

der Ausgangsbedingungen, 

- eine vierte Programmroutine zur rechnergestQtzten Anwendung des 
matrizentechnisch transformierten Minimierungskriteriums auf die Punkte der 
Begrenzungs-Punktmengen, 



f 6. Computer-Programm-Produkt beinhaltend einen maschinenlesbaren Pro- 
grammtrager, auf dem ein Computer-Programm nach Anspmch 5 in Form von 
elektronisch auslesbaren Steuersignalen gespeichert ist. 
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Zusammenfessung 

Verfahren zur rechnergestutzten Ermittlung einer treibstoffoptimalen Ansteueaing 
von Dusen 

Beschrieben wird ein Verfahren zur rechnergestutzten Ermittlung einer treibstoffopti- 
malen Ansteuerung von Dusen gemaB einer Ansteuervorschrift b = Ax. Dabei erfolgt 
rechnergestutzt eine definierte matrizentechnische Transformation von Ausgangsbe- 
dingungen fur den Massenfluss der DGsen und des Minimierungskriteriums, eine 
datenverarbeitungstechnische Darstellung einer gebmetrischen Beschreibung der 
matrizentechnisch transformierten Ausgangsbedingungen, eine rechnergestutzte 
Ermittlung von Begrenzungs-Punktmengen der geometrischen Beschreibung der 
Ausgangsbedingungen durch ein rechnergestutztes, geometrisches Suchverfahren 
im Vektorraum und die Anwendung des matrizentechnisch transformierten l\/linimie- 
rungskriteriums auf die Punkte der Begrenzungs-Punktmengen. 

(Fig. 2) 
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